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Description 

The invention relates to a method of determin- 
ing spectrum parameters of a spectrum in a fre- 
quency domain, related to spectroscopic signals in 
the time domain, in which method the spec- 
troscopic signals derived from a substance are 
sampled in order to obtain sampling values and are 
approximated by a complex model function which 
contains the spectrum parameters and exponen- 
tially damped sinusoids, in which method there is 
also made an initial estimate, based on the sam- 
pling values, for the spectrum parameters occurring 
non-linearly in the model function, the spectrum 
parameters being accurately estimated from the 
model function by iteration with a least-squares 
optimisation procedure, and prior knowledge is in- 
troduced into the model function. 

The invention also relates to a device for deter- 
mining spectrum parameters of a spectrum in a 
frequency domain, related to spectroscopic signals 
in the time domain, which device comprises 
means, for generating the spectroscopic signals in 
a substance, sampling means for obtaining sam- 
pling values from the spectroscopic signals, and 
display means for displaying the spectrum, and 
also comprises programmed means which com- 
prise a complex model function containing the 
spectrum parameters and exponentially damped 
sinusoids, and are suitable for storing prior knowl- 
edge of the spectroscopic signals as well as for 
making an initial estimate for the spectrum param- 
eters occurring non-linearly in the model function 
and for executing a least-squares optimisation pro- 
cedure for the model function on the basis of the 
sampling values, the programmed means further- 
more comprising iteration means for accurately es- 
timating, using iteration, the spectrum parameters 
by means of the least-squares optimisation proce- 
dure. 

A method of this kind is inter alia suitable for 
signals which are mainly exponentially damped, for 
example for signals obtained during a magnetic 
resonance experiment from an entire body as the 
substance as well as from a part of the body. The 
method can also be used, for example for X-ray 
spectroscopy or FT infrared spectroscopy. 

A method of this kind is disclosed in an article 
by J.W.C. van der Veen et al, "Accurate Quantifica- 
tion of in Vivo 31 P NMR Signals Using the Variable 
Projection Method and Prior Knowledge", Magnetic 
Resonance in Medicine, Vol. 6, No. 1, January 
1988, pp.92-98. 

The spectrum parameters of a spectrum re- 
lated to spectroscopic signals (signals obtained 
during an NMR experiment) are derived via the 
sampling values (represented in a time domain, 
which is in contrast with the spectrum which is 



represented in a frequency domain), directly in the 
time domain. A least-squares optimisation proce- 
dure is used for fitting the sampling values to a 
model function of spectroscopic signals, in which 

5 procedure prior knowledge concerning spectral 
components can be introduced. The method can in 
principle be used for an arbitrary model function. In 
said article, for example a model function is as- 
sumed which contains exponentially damped 

w sinusoids. By operation in the time domain, ar- 
bitrary parts of the spectroscopic signals can be 
omitted without giving rise to serious fitting prob- 
lems. For example, when in the case of NMR the 
substance contains immobile nuclei and mobile 

75 nuclei, the sampling values (of a quickly exponen- 
tially decaying spectroscopic signal) originating 
from the immobile nuclei can be simply omitted. 
The sampling values of the tail of an exponentially 
decaying spectroscopic signal can also be omitted. 

20 Thus, convolution effects associated with a trans- 
formation to the frequency domain in conjunction 
with weighting, such as base line (wide background 
in the spectrum due to quickly decaying signals 
from the immobile nuclei) and line-shape distor- 

25 tions (due to phase shifts) can be avoided. Fitting 
is performed partly non-iteratively in order to obtain 
an initial estimate for at least one of the spectrum 
parameters used as initial values for further itera- 
tive fitting. The non-iterative fitting is comparatively 

30 fast and the further iterative fitting enables the 
method to be used successfully also in the case of 
a poor signal-to-noise ratio. For parameters which 
occur linearly in the model function, such as com- 
plex amplitudes, no initial values are required, re- 

35 suiting in a faster method. Said article describes 
the method inter alia for a model function of ex- 
ponentially damped sinusoids in which parameters 
such as amplitude, damping factor and frequency 
occur. Using the least-squares optimisation proce- 

40 dure, the sampling values are iteratively fitted as 
well as possible to the model function, initial values 
being required for the damping factor and fre- 
quencies which occur as non-linear parameters in 
the model function. The prior knowledge may com- 

45 prise a time shift of an instant of a first sample with 
respect to a time origin and predetermined phases. 
The least-squares optimisation procedure can be 
written in the form of a matrix. This is represented 
by the formula 4 on page 94 of said article, show- 
so ing a least-squares solution for the amplitude. All 
terms in the inverse matrix occurring therein are 
successively calculated, like the terms in the prod- 
uct of the Hermitic matrix with the data vector. In 
the matrix form a large number of product terms 

55 occur, so that a substantial amount of computation 
time is required for executing the method by 
means of a computer; this aspect becomes more 
significant in the case of a comparatively large 
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number of sampling values. 

It is an object of the invention to provide a 
method which, when the model function contains 
exponentially damped sinusoids, requires substan- 
tially less computation time when executed by 
means of a computer. 

To achieve this, a method in accordance with 
the invention is characterized in that prior to the 
iteration the model function is linearised with re- 
spect to the spectrum parameters occurring non- 
linearly, after which inner product terms occurring 
in the least-squares procedure are converted into 
an analytic form which is used during iteration. The 
number of columns of the matrix is increased with 
respect to the matrix in said article due to residual 
parameters which remain after linearisation, for ex- 
ample frequency and damping. The invention is 
based on the idea to convert first the inner product 
terms into analytic form in order to solve the least- 
squares problem in the case of exponentially 
damped sinusoids as the model function, rather 
than calculate, like in the known method, all matrix 
products one term after the other for ail possible 
summations by means of row and column mul- 
tiplications. The computation time is thus substan- 
tially reduced. In comparison with the known meth- 
od, the method in accordance with the invention is 
30 times faster in the case of 17 sinusoids and 512 
sampling values. Random tests have shown that 
the Cramer-Rao lower limit for estimating the pa- 
rameters is at least substantially reached, so that 
the method produces reliable parameters. The 
computation time increases in proportion to the first 
power of the number of sampling values, assuming 
that the other parameters, for example the number 
of sinusoids, remain the same. 

A version of a method in accordance with the 
invention is characterized in that terms of the ana- 
lytic form are recursively determined. Thus, the 
method can be simply implemented in a computer. 

A version of a method in accordance with the 
invention is characterized in that the model function 
contains a polynomial whereby the exponentially 
damped sinusoids are multiplied. A sum of ex- 
ponentially damped sinusoids does not always con- 
stitute a good model function in all circumstances. 
By multiplication by a polynomial, a model function 
is obtained which is adapted to circumstances. The 
inner product terms retain their comparatively sim- 
ple form. 

A version of a method in accordance with the 
invention is characterized in that the initial estimate 
of the spectrum parameters occurring non-linearly 
in the model function is made by means of a non- 
iterative estimation procedure. Initial values can 
thus be determined for parameters which occur 
non-linearly in the model function. For a non-itera- 
tive estimation procedure HSVD (Hankel matrix sin- 



gular value decomposition) can be used as de- 
scribed in "Improved Algorithm for Noniterative 
Time-Domain Model Fitting to Exponentially 
Damped Magnetic Resonance Signals", H. Bark- 

5 huijsen et al, Journal of Magnetic Resonance 73, 
pp. 553-557, 1987. 

A preferred version of a method in accordance 
with the invention is characterized in that the initial 
estimate of the spectrum parameters occurring 

w non-linearly in the model function is made on the 
basis of a coarse spectrum obtained from the sam- 
pling values by Fourier transformation. In the case 
of a comparatively poor signal-to-noise ratio, first a 
filtering operation can be performed. Initial values 

75 of the parameters which occur non-linearly in the 
model function can be determined from the coarse 
spectrum. These parameters, frequency and damp- 
ing, can be determined directly from a display 
screen on which the coarse spectrum is displayed 

20 by picking up, using a so-called mouse, at peak 
locations in the coarse spectrum the peak locations 
themselves (frequencies) and the peak width at half 
peak height (a measure of the damping) for supply 
to the programmed means in a computer. When 

25 the coarse spectrum contains excessive phase dis- 
tortion, first a phase correction method incorpo- 
rated in the programmed means can be executed, 
possibly after filtering (in the case of a compara- 
tively poor signal-to-noise ratio), for example as 

30 described in an article by C.H.Sotak et al in Journal 
of Magnetic Resonance, Vol. 57, pp. 453-462, 
1984. The phase correction of the coarse spectrum 
and the determination of the initial values can also 
be performed by means of a method as described 

35 in EP-A 0 316 041 , which forms part of the prior art 
under Art. 54(3) EPC, "Method and device for 
automatic phase correction of complex NMR spec- 
tra", by J. van Vaals. 

The invention will be described in detail 

40 hereinafter with reference to a drawing, therein 

Fig. 1 diagram matically shows a device in ac- 
cordance with the invention, 
Fig. 2 shows, in the form of a matrix, a set of 
linear equations in accordance with the inven- 

45 tion, to be used for each iteration step during 
the iteration process, 

Fig. 3 shows the analytic form to be substituted, 
in accordance with the invention, in the set of 
linear equations given in Fig. 2, 
so Fig. 4 shows a spectrum obtained by means of 
the method in accordance with the invention, 
and 

Fig. 5 shows a table containing parameters of 

the spectrum of Fig. 4. 
55 Fig. 1 diagrammatically shows a device 1 in 

accordance with the invention. Specifically, a mag- 
netic resonance device for obtaining spectroscopic 
signals is shown by way of example. The device 1 



5 



EP 0 341 783 B1 



6 



comprises magnet coils 3 for generating a steady, 
uniform magnetic field Bo which are arranged with- 
in a shielded space 2. Within the magnet coils 3 
there may be arranged a substance. The substance 
is, for example a body containing solid and liquid 
components. During operation of the device, the 
substance being arranged within the magnet coils 
3, a slight excess of nuclear spins (of nuclei having 
a magnetic moment) will be oriented in the same 
direction as the field B 0 . From a macroscopic point 
of view this is to be considered as a magnetisation 
M, i.e. an equilibrium magnetisation. The magnet 
coils 3 may be resistance magnets or supercon- 
ducting magnets; the device then comprises a d.c. 
power supply source 4 for powering the magnet 
coils 3. When the magnet coils 3 are constructed 
as permanent magnets, the d.c. power supply 
source 4 will be omitted. The device 1 furthermore 
comprises means 5 for generating spectroscopic 
signals, which means include a reference generator 
6 for generating a reference signal for supply to an 
r.f. modulator 7. When activated, the r.f. modulator 
7, supplied an r.f. electromagnetic pulse to an r.f. 
power amplifier 8 which amplifies the r.f. elec- 
tromagnetic pulse for supply to a transmitter coil 9 
which is arranged near the magnet coils 3. When 
the substance arranged within the magnet coils is 
irradiated, the equilibrium of the magnetisation M is 
disturbed and after irradiation there will be gen- 
erated a magnetic resonance signal when the mag- 
netisation returns to the state of equilibrium, a so- 
called FID signal (Free Induction Decay). The reso- 
nance signal contains spectroscopic information; 
the signal will be referred to hereinafter as a spec- 
troscopic signal. The device 1 also comprises 
means 10 for receiving the spectroscopic signals, 
including a receiver coil 11 for receiving the spec- 
troscopic signals which is connected to a detector 

12 for the detection of the modulated spectroscopic 
signal, for example by quadrature detection. The 
detector 12, being coupled to the reference gener- 
ator 6, comprises low-pass filters and analog-to- 
digital converters as the sampling means for sam- 
pling the spectroscopic signals received and de- 
tected. The device 1 also comprises control means 

13 which serve inter alia for control and timing of 
the r.f. generator 7. The device furthermore com- 
prises processing means 14 for processing digitis- 
ed spectroscopic signals 15 (sampling values). The 
processing means 14 are coupled to the control 
means 13. The processing means 14 comprise 
programmed means 16 inter alia for making an 
initial estimate for at least one spectrum parameter 
of a spectrum related to the spectroscopic signals 
in order to perform a least-squares optimisation 
procedure on the basis of the sampling values. The 
programmed means 16 furthermore contain a 
model function of exponentially damped sinusoids, 



a non-iterative estimation procedure, iteration 
means and Fourier transformation means (not 
shown in detail). The processing means 14 also 
comprise a RAM memory 17 which is coupled to 

5 the programmed means 16. The device 1 further- 
more comprises display means 18, inter alia for the 
display of a spectrum of the spectroscopic signals. 
A cursor control device 19 is coupled to the display 
means 18 and the processing means 14. The cur- 

w sor control device may be a so-called mouse or 
track ball; as is known, the mouse is formed by a 
carriage provided with a motion sensor which out- 
put signals in response to motion of the mouse. 
The signals of the motion sensor are applied to the 

75 processing means 14 which move a cursor as an 
indicator across a display screen of the display 
device in dependence of the signals of the motion 
sensor. When a given screen position is reached, 
the screen position can be applied to the process- 

20 ing means by operation of a control member on the 
mouse. Inter alia for spectroscopy of a part of the 
substance, so-called volume-selective spectro- 
scopy, the device 1 comprises gradient coils 20 
which are controlled by a gradient coil control 

25 device 21. The gradient coil control device 21 is 
controlled by the control means 13. The gradient 
coils 20 serve to generate field gradients which are 
superposed on the field B 0 . Generally, there can be 
generated three gradients Gx, Gy and Gz whose 

30 field direction coincides with the direction of the 
steady, uniform magnetic field Bo and the respec- 
tive gradient directions of which extend mutually 
perpendicularly. For a more detailed description of 
a magnetic resonance device, reference is made to 

35 the handbook "Practical NMR Imaging", by M.A. 
Foster and J. M.S. Hutchinson, IRL Press, 1987, 
ISBN 1-85221-011-7, Chapter 1, pages 1-48 which 
describes the general principles of NMR and 
shows an NMR device. In the case of spectroscopy 

40 of the substance or part of the substance, the aim 
is to obtain information from the substance in the 
form of a spectrum (in the frequency domain). For 
NMR, where nuclear spins are excited in the 
steady field B 0 , the relation w =gamma.(1-a).B 0 

45 holds for excitation of nuclear spins, where o is a 
resonant frequency of the nuclear spins, gamma is 
a gyromagnetic ratio, and a is a so-called shielding 
constant of a nucleus. Spin resonance will occur 
when the resonant frequency is applied. From a 

50 macroscopic point of view this implies, vis-a-vis a 
stationary observer, a precessional motion of the 
magnetisation M around the field B 0 at the reso- 
nant frequency. When the substance contains 
nuclei having a different a, the resonant frequency 

55 will also be different for these nuclei. This effect is 
referred to as chemical shift and becomes manifest 
as different resonant peaks in a spectrum obtained 
from the resonance signal. In the time domain the 
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resonance signal can in principle be described as a 
superposition of exponentially damped sinusoids. 
The chemical shifts are due to the fact that nuclei 
effectively experience a different magnetic field 
due to a different shielding under the influence of 
the presence of neighbouring electrons. The neigh- 
bouring electrons shield the nuclei to a given de- 
gree from the field B 0 . Thus, chemical components 
can be identified in a substance by means of NMR 
spectroscopy. Depending on a selected frequency 
of the r.f. magnetic pulse, a given spectrum will be 
obtained. For example, a phosphor spectrum ( 31 P), 
a carbon spectrum ( 13 C) or a proton spectrum ( 1 H) 
can be obtained. The surface area below the reso- 
nance peaks in the spectrum is directly propor- 
tional to the number of nuclei making a contribution 
for a given frequency. The surface area corre- 
sponds to the amplitude in the time domain. In 
addition to chemical shifts, another interaction also 
has an effect on the spectrum, the so-called spin- 
spin coupling, which causes the splitting off of a 
resonance line which is usually small in compari- 
son with the chemical shift. For a more detailed 
description of NMR spectroscopy reference is 
made to pages 249-274 of said handbook by Fos- 
ter and Hutchinson. 

Fig. 2 shows a set of linear equations (2A) to 
be used in accordance with the invention, for each 
iteration step during the iteration process. After an 
initial estimate of at least one spectrum parameter, 
frequency and damping of peaks occurring in the 
spectrum in the present example, the parameters 
frequency and damping and also amplitude are 
accurately estimated in a number of iteration steps 
on the basis of the set (2A). Therein, the column 
vector c relates to the amplitudes of the peaks of 
the spectrum and Re c' represents the residual 
parameters Ap of the damping p and Im c' repre- 
sents the residual parameters Anu of the frequency 
nu, M1 to M6 being matrices in the partitioned form 
of the extended Fisher data matrix M of the least- 
squares problem, extended with columns for the 
residual parameters A p and Anu for all peaks, and 
vi to V3 a representation of the sampling values. 
The set (2A) is formed by linearisation of the least- 
squares problem. The procedure is based on a 
model function for the magnetic resonance signal 
x n of the form: 

x n = I k c k .exp[( Pk + 27ri.nu k )(n+a) + i^ k ] (1) 

where x n is a model function of the signal x n , c k , 
nu k and 

p k are parameters of the spectrum of x n , being 
amplitude, 

frequency and damping, respectively, for k=1 

K, in which K is the number of peaks in the model 
function. Furthermore, i = (-1) 1/2 . First an initial es- 



timate is made for p k and nu k . In the working point 
(po. nuo) obtained for the parameters p and nu, for 
each k x n is linearised, producing terms of the 
form: 

5 

[c + (n + 5)c'].exp[(p 0 + 2?ri.nuo)(n + 5) + i<£o ] (2) 

in x n .p= po+Ap and nu=nuo+Anu. Therein, n 
represents a 

10 sampling time of the resonance signal, n=0, 1, .... 
N-1 for N sampling operations at instants t n , 5 
denotes a time shift caused by instrumental limita- 
tions and/or the deliberate omission of sampling 
values, </> k denotes a phase shift. 5 and £ k can be 

15 introduced into the model x n as prior knowledge. 
Furthermore, after linearisation Re(c')=cAp and 
lm(c') = 2 TrcAnu, which can be simply dem- 
onstrated. Formulas (2B) to (2G) follow from further 
elaboration of the least-squares problem after 

20 linearisation. In the formule (2A) to (2G), " is a 
transposition and * is a complex conjugated opera- 
tion. f k (t n ) is the k th model function. The information 
matrix is dimensioned to 3Kx3K. Furthermore: 

25 f k (t n ) = exp[( Pk + 27ri.nu k )(n + 5) + i* k ] (3) 

After each iteration step the current values of p k 
and nu k are replaced, utilising the current values of 
c k . When c' is sufficiently small, the iteration stops. 

30 A spectrum is obtained by transforming the model 
obtained after iteration to the frequency domain, for 
example by means of the Fourier transformation 
means. The spectrum can also be calculated di- 
rectly from the spectrum parameters obtained, us- 

35 ing known Lorentz formules, so that any artefacts 
due to transformation are avoided. The spectrum is 
displayed by means of the display means 18. The 
spectrum parameters can also be represented in 
the form of a table; in that case the Fourier trans- 

40 formation may be omitted. Random tests have 
demonstrated that the so-called Cramer-Rao lower 
limit is at least substantially reached. For a further 
description of the Cramer-Rao lower limit in relation 
to the Fisher data matrix reference is made to 

45 "Handbook of Measurement Science", Vol. 1, P.H. 
Sydenham, 1982, Wiley, pp. 335-339. For a de- 
scription of the generally known Fourier transforma- 
tion means reference is made to the handbook "An 
Introduction to the Analysis and Processing of Sig- 

50 nals" P.A. Lynn, McMillan 1983, pp. 57-64. 

Fig. 3 shows the analytic form for the inner 
product terms to be substituted in the set of linear 
equations shown in Fig. 2. By elaborating in the 
formulas (2B) to (2G) the product f k (t n )V(t n ), multi- 

55 plied by a power of t n , for f k (t n ) in accordance with 
(3), terms according to (3A) in Fig. 3 are produced 
in as far as the part thereof which is dependent on 
the sampling instant is concerned; therein, 1 in- 



9 



EP 0 341 783 B1 



10 



dicates the formulas corresponding to the formulas 
(2B), (2C) and (2D) in Fig. 2. In formula (3A): 

z = exp[p k + Pk < + 27ri(nu k -nu k ')] (4) 

Formulas (3B) and (3C) show recursive relations for 
the respective terms corresponding to l = 0, 1 and 
2. When z approximates 1, an inaccuracy may 
occur when use is made of the recursive formules 
(3B) and (3C). By introducing the variable eps = z- 
1, this problem an be circumvented. (3D) to (3F) 
show formules for use in such a case. The com- 
putation time required will then slightly increase 
again, notably in the case of a large N, because of 
the combinations occurring therein, for example 
when N exceeds n + 1 . Without substantial loss of 
accuracy truncation can take place for a suitable 
value of n, so that the computation efficiency is 
increased again. For example, in the case of 512 
sampling values it will suffice to truncate the sum- 
ming operations for, for example, n = 10. 

Fig. 4 shows a spectrum obtained by means of 
the method in accordance with the invention for in 
vivo 13 C NMR spectroscopy of the liver of a rat. 
The number of complex sampling operations 
amounted to 512 with a sampling interval of 50 us. 
The frequency scale is given in kHz relative to the 
irradiation frequency (the chemical shift appears as 
a spectrum). In Fig. 4 the reference 4A is the 
cosine FFT (Fast Fourier Transform) of the sam- 
pling values ( a window is used in order to avoid 
artefacts), the reference 4B is a spectrum obtained 
by means of a method other than the method in 
accordance with the invention, and the reference 
4C is the cosine FFT of the model which has been 
fitted in the time domain by means of the method 
in accordance with the invention and in which the 
prior knowledge concerning the phases has been 
introduced. In accordance with the invention, one 
operation step required 15 seconds, using the 
same device, whereas the known method required 
a period of time which is 30 times longer. The 
references p1 tp p17 denote resonance peaks. The 
spectrum shown in Fig. 4C was obtained by means 
of 8 iteration steps. 

Fig. 5 shows a table containing the parameters 
of the spectrum 4c in Fig. 4. The number of 
exponentially damped sinusoids was initially set at 
25; it was found that 17 sinusoids were of spec- 
troscopic relevance. For each resonance peak 
there are shown parameters obtained using a 
method other than the method in accordance with 
the invention, in this case by means of HSVD 
without utilising prior knowledge and therebelow 
the parameters obtained utilising the method in 
accordance with the invention with introduction of 
prior knowledge. The errors indicated amount to 
twice the standard deviation. 



Many alternatives are feasible for those skilled 
in the art without departing from the scope of the 
invention. For example, in the case of NMR many 
pulse sequences can be used for spectroscopy in 

5 order to obtain resonance signals, both volume 
selective as well as non-volume selective. For a 
volume-selective pulse sequence use can be made 
of the SPARS method, as disclosed in Magnetic 
Resonance Imaging, Vol. 4, pp.237-239, 1986. 

70 Even more prior knowledge can be used, for exam- 
ple information concerning multiplets in the spec- 
trum. For a given signal-to-noise ratio for in vivo 
NMR, the use of prior knowledge will often be 
required in order to achieve suitable quantification. 

15 

Claims 

1. A method of determining spectrum parameters 
of a spectrum in a frequency domain related to 

20 spectroscopic signals in the time domain, in 

which method the spectroscopic signals de- 
rived from a substance are sampled in order to 
obtain sampling values and are approximated 
by a complex model function which contains 

25 the spectrum parameters and exponentially 

damped sinusoids, in which method there is 
also made an initial estimate, based on the 
sampling values, for the spectrum parameters 
occurring non-linearly in the model function, 

30 the spectrum parameters being accurately es- 

timated from the model function by iteration 
with a least-squares optimisation procedure, 
and prior knowledge is introduced into the 
model function, characterized in that prior to 

35 the iteration the model function is linearised 

with respect to the spectrum parameters oc- 
curring non-linearly, after which inner product 
terms occurring in the least-squares procedure 
are converted into an analytic form which is 

40 used during iteration. 

2. A method as claimed in Claim 1, characterized 
in that terms of the analytic form are recur- 
sively determined. 

45 

3. A method as claimed in Claim 1 or 2, char- 
acterized in that the model function contains a 
polynomial whereby the exponentially damped 
sinusoids are multiplied. 

50 

4. A method as claimed in Claim 1, 2 or 3, 
characterized in that the initial estimate of the 
spectrum parameters occurring non-linearly in 
the model function is made by means of a 

55 non-iterative estimation procedure. 

5. A method as claimed in Claim 1, 2 or 3, 
characterized in that the initial estimate of the 
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spectrum parameters occurring non-linearly in 
the model function is made on the basis of a 
coarse spectrum obtained from the sampling 
values by Fourier transformation. 

6. A device for determining spectrum parameters 
of a spectrum in a frequency domain related to 
spectroscopic signals in the time domain, 
which device comprises means (3, 5), for gen- 
erating the spectroscopic signals in a sub- w 
stance, sampling means (12) for obtaining 
sampling values from the spectroscopic sig- 
nals, and display means (18) for displaying the 
spectrum, and also comprises programmed 
means (16) which comprise a complex model 75 
function containing the spectrum parameters 

and exponentially damped sinusoids, and are 
suitable for storing prior knowledge of the 
spectroscopic signals as well as for making an 
initial estimate for the spectrum parameters 20 
occurring non-linearly in the model function 
and for executing a least-squares optimisation 
procedure for the model function on the basis 
of the sampling values, the programmed 
means (16) furthermore comprising iteration 25 
means for accurately estimating, using iter- 
ation, the spectrum parameters by means of 
the least-squares optimisation procedure, char- 
acterized in that the iteration means contain 
the model function linearised with respect to 30 
the spectrum parameters occurring non-linear- 
ly, inner product terms occurring in the least- 
squares optimisation procedure being present 
in an analytic form, which analytic form serves 
for use in the iteration means. 35 

7. A device as claimed in Claim 6, characterized 
in that the programmed means (16) also con- 
tain a non-iterative estimation procedure for 
making an initial estimate of the spectrum pa- 40 
rameters occurring non-linearly in the model 
function. 

8. A device as claimed in Claim 6 which com- 
prises Fourier transformation means, character- 45 
ized in that the programmed means (16) are 

also suitable for determining, in conjunction 
with the Fourier transformation means, a 
coarse spectrum from the sampling values and 
for displaying the coarse spectrum by way of 50 
the display means (18), the device comprising 
a cursor control device (19) which is coupled 
to the display device (18) in order to pick up 
parameters from the coarse spectrum as an 
initial estimate for the spectrum parameters 55 
occurring non-linearly in the model function. 
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Patentanspruche 

1. Verfahren zum Bestimmen von Spektralpara- 
metern eines auf spektralanalytische Signale in 
der Zeitebene bezogenen Spektrums in einer 
Frequenzebene, bei dem die von einer Sub- 
stanz abgeleiteten spektralanalytischen Signale 
zum Erhalten von Abtastwerten abgetastet und 
durch eine komplexe Modellfunktion angena- 
hert werden, die die Spektral parameter und 
exponentiell gedampfte Sinuskurven enthalt, 
und bei dem auflerdem eine anfangliche 
Schatzung fur die nicht-linear in der Modell- 
funktion auftretenden Spektralparameter auf 
der Basis der Abtastwerte vorgenommen wird, 
wobei die Spektralparameter, von der Modell- 
funktion ausgehend, durch Iteration mit einer 
Optimierung nach der Methode der kleinsten 
Quadrate sehr genau geschatzt werden und 
Vorkenntnisse in die Modellfunktion eingefuhrt 
werden, dadurch gekennzeichnet, daB die Mo- 
dellfunktion vor der Iteration in bezug auf die 
nicht-linear auftretenden Spektralparameter li- 
nearisiert wird und anschlieSend in der Metho- 
de der kleinsten Quadrate auftretende 
Skalarprodukt-Terme in eine analytische Form 
umgewandelt werden, die wahrend der Itera- 
tion verwendet wird. 

2. Verfahren nach Anspruch 1, dadurch gekenn- 
zeichnet, daB Terme der analytischen Form 
rekursiv bestimmt werden. 

3. Verfahren nach Anspruch 1 oder 2, dadurch 
gekennzeichnet, daB die Modellfunktion ein Po- 
lynom enthalt, mit dem die exponentiell ge- 
dampften Sinuskurven multipliziert werden. 

4. Verfahren nach Anspruch 1, 2 oder 3, dadurch 
gekennzeichnet, daB die anfangliche Schat- 
zung der nicht-linear in der Modellfunktion auf- 
tretenden Spektralparameter mit Hilfe einer 
nicht-iterativen Schatzmethode erfolgt. 

5. Verfahren nach Anspruch 1 , 2 oder 3, dadurch 
gekennzeichnet, daB die anfangliche Schat- 
zung der nicht-linear in der Modellfunktion auf- 
tretenden Spektralparameter auf der Basis ei- 
nes groben Spektrums durchgefuhrt wird, das 
durch Fourier-Transformation aus den Abtast- 
werten erhalten worden ist. 

6. Anordnung zum Bestimmen von Spektralpara- 
metern eines auf spektralanalytische Signale in 
der Zeitebene bezogenen Spektrums in einer 
Frequenzebene, wobei die Anordnung Mittel 
(3, 5) zum Erzeugen der spektralanalytischen 
Signale in einer Substanz, Abtastmittel (12) 
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zum Erhalten von Abtastwerten aus den spek- 
tralanalytischen Signalen, und Anzeigemittel 
(18) zur Darsteliung des Spektrums umfaBt und 
auch programmierte Mittel (16) enthalt, die 
eine komplexe Modellfunktion mit den Spek- 
tralparametern und exponentiell gedampften 
Sinuskurven umfassen und sich zum Spej- 
chern der Vorkenntnisse uber die spektralana- 
lytischen Signale eignen sowie zum Durchfuh- 
ren einer anfanglichen Schatzung fUr die nicht- 
linearen Spektralparameter in der Modellfunk- 
tion und zum Durchfuhren einer Optimierung 
nach der Methode der kleinsten Quadrate fur 
die Modellfunktion auf der Basis der abgetaste- 
ten Werte, wobei die programmierten Mittel 
(16) auBerdem Iterationsmittel zum genauen 
Schatzen der Spektralparameter durch Iteration 
mit Optimierung nach der Methode der klein- 
sten Quadrate enthalt, dadurch gekennzeich- 
net, daB die Iterations-Mittel die in bezug auf 
die nicht-linear auftretenden Spektralparameter 
linearisierte Modellfunktion enthalten, wobei 
Skalarprodukt-Terme, die bei der Optimierung 
nach der Methode der kleinsten Quadrate auf- 
treten, in einer analytischen Form vorliegen, 
die fur die Verwendung in den Iterations-Mit- 
teln dient. 

7. Anordnung nach Anspruch 6, dadurch gekenn- 
zeichnet, dafl die programmierten Mittel (16) 
auch eine nicht-iterative Schatzungsprozedur 
fOr die DurchfOhrung einer anfanglichen Schat- 
zung fur die in der Modellfunktion nicht-linear 
auftretenden Spektralparameter enthalten. 

8. Anordnung nach Anspruch 6 mit Fourier- 
Transformations-Mitteln, dadurch gekennzeich- 
net, daB die programmierten Mittel (16) sich in 
Verbindung mit den Fourier-Transformations- 
Mitteln auch fur die Bestimmung eines groben 
Spektrums aus den Abtastwerten und zur An- 
zeige des groben Spektrums mit Hilfe der An- 
zeigemittel (18) eignen, wobei die Anordnung 
eine Cursor-Steuereinrichtung (19) enthalt, die 
mit der Anzeigeeinrichtung (18) gekoppelt ist, 
urn die Parameter aus dem groben Spektrum 
als anfangliche Schatzung fUr die nicht-linear 
in der Modellfunktion auftretenden Spektralpa- 
rameter aufzunehmen. 

Revendications 

1. Procede peur determiner les parametres spec- 
traux d'un spectre (1) dans un domaine de 
frequences en relation avec des signaux spec- 
troscopiques (2) dans le domaine temporel, 
procede dans lequel les signaux spectroscopi- 
ques derives d'une substance sont echantillon- 



nes aux fins d'obtenir des valeurs d'echantil- 
lonnage et sont estimes approximativement 
par une fonction modele complexe qui contient 
les parametres spectraux et des sinusoVdes 

5 amorties exponentiellement, procede dans le- 

quel est egalement faite une estimation initiale, 
bas6e sur les valeurs d'^chantillonnage, pour 
les parametres spectraux apparaissant non li- 
n£airement dans la fonction modele, les para- 

w metres spectraux etant estimes avec precision 

d'apres la fonction modele par iteration au 
moyen d'une procedure d'optimalisation par 
les moindres carres, et la connaissance ante- 
rieure est introduite dans la fonction modele, 

15 caracterise en ce qu'avant I'iteration, la fonc- 

tion modele est linearisee par rapport aux pa- 
rametres spectraux apparaissant non lineaire- 
ment, apres quoi les termes, qui sont des 
produits internes apparaissant dans I'operation 

20 par les moindres carres, sont convertis en une 

forme analytique qui est utilisee pendant I'ite- 
ration. 

2. ProcedS suivant la revendication 1, caracterise 
25 en ce que les termes de la somme analytique 

sont determines par recurrence. 

3. Procede suivant la revendication 1 ou 2, carac- 
terise en ce que la fonction modele contient un 

30 polynome au moyen duquel les sinusoVdes 

amorties exponentiellement sont multiplies. 

4. Procede suivant la revendication 1, 2 ou 3, 
caracterise en ce que Testimation initiale des 

35 parametres spectraux apparaissant non lineai- 

rement dans la fonction modele est effectuee 
au moyen d'une procedure d'estimation non 
interative. 

40 5. Procede suivant la revendication 1, 2 ou 3, 
caracterise en ce que Testimation initiale des 
parametres spectraux apparaissant non lineai- 
rement dans la fonction modele est effectuee 
sur la base d'un spectre grossier obtenu a 

45 partir des valeurs d'echantillonnage par trans- 

formation de Fourier. 

6. Dispositif pour determiner les parametres 
spectraux d'un spectre (1) dans un domaine 

so de frequences en relation avec des signaux 

spectroscopiques (2) dans le domaine tempo- 
rel, lequel dispositif comprend des moyens (3, 
5) pour engendrer les signaux spectroscopi- 
ques dans une substance, un moyen d'£chan- 

55 tillonnage (12) pour obtenir des valeurs 

d'echantillonnage a partir des signaux spec- 
troscopiques et un moyen d'affichage (18) 
pour afficher le spectre, et comprend aussi un 
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moyen programme (16) qui comprend une 
fonction module complexe contenant les para- 
metres spectraux et des sinusoTdes amorties 
exponentiellement et qui convient pour mettre 
en memoire la connaissance anterieure des 5 
signaux spectroscopiques, de meme que pour 
faire une estimation initiale des param&tres 
spectraux apparaissant non lineairement dans 
la fonction modele et pour ex£cuter une opera- 
tion d'optimalisation par les moindres carres 10 
pour la fonction modele sur la base des va- 
leurs d'£chantillonnage, le moyen programme 
(16) comprenant, de surcroit, un moyen d'ite- 
ration pour estimer avec precision, par itera- 
tion, les parametres spectraux a I'aide de la 75 
procedure d'optimalisation par les moindres 
carres, caract£rise en que le moyen d'iteration 
contient la fonction modele linearisee par rap- 
port aux parametres spectraux apparaissant 
non lineairement, les termes qui sont des pro- 20 
duits internes apparaissant dans la procedure 
d'optimalisation par les moindres carres etant 
presents sous une forme analytique, laquelle 
forme analytique sert dans le moyen dera- 
tion. 25 

7. Dispositif suivant la revendication 6, caract£ri- 
se en ce que le moyen programme (16) 
contient aussi une procedure d'estimation non 
iterative pour effectuer une estimation initiale 30 
des parametres spectraux apparaissant lineai- 
rement dans la fonction module. 

8. Dispositif suivant la revendication 6, qui com- 
prend un moyen de transformation de Fourier, 35 
caracterise en ce que le moyen programme 

(16) convient aussi pour determiner, conjointe- 
ment avec le moyen pour transformation de 
Fourier, un spectre grossier h partir des va- 
leurs d'£chantillonnage et pour afficher le 40 
spectre grossier a I'aide du moyen d'affichage 
(18), le dispositif comprenant a dispositif de 
commande de curseur (19) qui est couple au 
moyen d'affiche (18) afin de collector les para- 
metres du spectre grossier en tant qu'estima- 45 
tion initiale des parametres spectraux apparais- 
sant non lineairement dans la fonction modele. 
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